Are stress-free membranes really 'tensionless' ? 
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In recent years it has been argued that the tension parameter driving the fluctuations of fluid 
membranes, differs from the imposed lateral stress, the 'frame tension'. In particular, stress-free 
membranes were predicted to have a residual fluctuation tension. In the present paper, this argument 
is reconsidered and shown to be inherently inconsistent - in the sense that a linearized theory, 
the Monge model, is used to predict a nonlinear effect. Furthermore, numerical simulations of 
one-dimensional stiff membranes are presented which clearly demonstrate, first, that the internal 
'intrinsic' stress in membranes indeed differs from the frame tension as conjectured, but second, 
that the fluctuations are nevertheless driven by the frame tension. With this assumption, the 
predictions of the Monge model agree excellently with the simulation data for stiffness and tension 
values spanning several orders of magnitude. 

PACS numbers: 87.16.dj, 68.03.Kn, 68.35.Md 



Fluid membranes are indispensable constituents of all 
living things Q . They are made of self- assembling am- 
phiphilic molecules, mostly lipids, which aggregate to 
bilayer sheets in an aqueous environment. The funda- 
mental properties of self-assembled membranes have been 
studied for many decades by experimental studies of sim- 
ple model membranes, by random interface theories, and 
by simulations on length scales ranging from the atomic 
to the mesoscopic scale. Whereas in nature and in typi- 
cal experimental setups (vesicles, supported membranes), 
membranes are usually under (slight) tension, most the- 
oretical and simulation studies have focused on the ideal 
'tensionless' case. 

In the present paper we address a basic question: What 
is the meaning of 'tensionless' in a membrane? Do 
theories and simulations/experiments really consider the 
same object if they refer to tensionless membranes? This 
has recently been questioned by a number of authors 0- 

@. 

Experimentally and in molecular simulations, the 
membrane tension is identified with the lateral stress, 
which is the mechanically accessible quantity. It corre- 
sponds to the force which the membrane would exert on a 
surrounding frame, hence it is also called 'frame tension' 
ct,. Membranes are said to be 'tensionless' if er f = 0. In 
fluctuating membranes, the frame tension is not neces- 
sarily identical with the internal membrane tension <r int , 
i.e., the conjugate variable to the total area in the free 
energy (at fixed number of lipids) , since the local mem- 
brane normal fluctuates. 

Theoretically, fluctuating membranes are commonly 
described in terms of effective interface models: The 
membranes are represented by elastic sheets, subject to 
an effective potential that depends on their conforma- 
tion Q . The simplest and most famous example of such 
a model is the Hclfrich Hamiltonian |8f 



for fluid membranes. Here the integral A A runs over the 
surface of the sheets, H and G are the total and the 
Gaussian curvature (i.e., the sum and the product of 
the inverse principal curvature radii), and k and k are 
the bending rigidity and the Gaussian rigidity, respec- 
tively. The first term with the 'bare tension' a couples 
directly to the total surface area A of the membranes. 
For fluctuating interfaces (e.g., between two coexisting 
fluid phases), this term dominates the effective interface 
Hamiltonian In theoretical studies of 'tensionless' 
membranes, cr is often set to zero. The bare tension 
is related, but again not necessarily identical with the 
"fluctuation tension" cr fluc , i.e., the coefficient of q 2 in 
the fluctuation spectrum of a planar membrane floj|. 

The relation between the externally imposed 'frame 
tension' a { and the internal tension a iat has intrigued sci- 
entists for quite some time. In this context, it is im- 
portant to note that Eq. JT]) leaves room for different 
interpretations. If ([TJ is taken to describe strictly incom- 
pressible membranes with fixed lipid but fluctuating 
number of lipids, a Q is proportional to the chemical po- 
tential jUo for adding a lipid [ll| . The number of degrees 
of freedom N then fluctuates together with the total area 
A (at fixed A/N). This is the point of view taken by Cai 



dA{aa + 2_H 2 + ^G} 
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et al. in a seminal work of 1994 At fixed projected 
area A p (but fluctuating number of lipids), they argued 
that the fluctuation tension cr fluc rcnormalizes to a anc = a t 
in the thermodynamic limit. More recent authors 0-@| 
have used a Hamiltonian of the form ([1]) to describe com- 
pressible membranes with fixed number of lipids. In that 
case, the number of degrees of freedom N is constant 
and <7 is the internal tensionQ c int . This situation has 
first been considered by Farago and Pincus in 2003/04, 
who concluded again tr fluc = a { for membranes with fixed 
projected area US]. 

This conclusion was recently questioned by Imparato 
[1] and Fournier and Barbetta |f3|. Based on free energy 
calculations as well as direct calculations of the thermally 
averaged lateral stress tensor [j| for almost planar inter- 
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faces (Monge representation, see below), they recover an 
earlier result by Farago and Pincus regarding the differ- 
ence between the internal tension er int and the frame ten- 
sion (TfQ, but they argue that the fluctuations are driven 
by cr int . The fluctuation tension is predicted to be always 
larger than the frame tension. Most notably, stress-free 
membranes with o~ { = are predicted to have a residual 
fluctuation tension er fluc (°) = k B TN/A, where N is the 
number of fluctuation degrees of freedom. The residual 
tension can be attributed to the contribution of the fluc- 
tuations on the interfacial free energy. On large length 
scales, it should dominate the fluctuation spectrum. 

Experimentally, however, the evidence for such an ef- 
fect is scarce. In experiments, the effect of fluctuations 
can be measured by monitoring the relative excess area 
((A — Ap)/A p ), and the frame tension a t can be imposed, 
e.g., by the Laplace pressure in micropipette experiments 
[12L 1 1 31 ] . The experimental findings in such setups were 
consistent with the assumption er fluc = a t . A difference 
between cr fluc and a { was only reported in one recent study 
of giant vesicles jTij], where er ftuc was inferred from the 
amplitude of fluctuations and <7 f from the shape of the 
vesicles on surfaces. 

Likewise, most simulations of fluctuating stress-free 
membranes [l5 - 21 1 or membrane stacks 22 1 could be de- 
scribed very satisfactorily by elastic theories for er auc = 0, 
i.e., the long- wavelength behavior is apparently domi- 
nated by ben ding modes. This holds for simulations of 
atomistic la, ll7| as well as coarse-grained la, Il8l422 
models, with two exceptions: In simulations of mem- 
branes under varying stress, Imparato [3| and Stecki 23 1 
independently report the fluctuation tension cr nuc to be 
larger than the frame tension a s over a range of frame 
tensions, in agreement with the theoretical prediction. In 
contrast, a more recent coarse-grained simulation study 
by Neder et al. dij yielded cr nuc = er t within the error at 
moderate frame tensions and <7 fluc < a f at extreme ten- 
sions. A reduction of cr fluc compared to a { cannot be ex- 
plained by any available theory. On the other hand, the 
membrane was also found to exhibit substantial struc- 
tural changes in this high tension regime, which presum- 
ably accounts for part of the effect. 

In sum, the majority of experimental and numerical 
evidence so far points towards er fluc = er f at low tensions, 
and against the existence of a residual fluctuation tension 
in stress- free membranes. On the other hand, neither 
the simulation data nor the experimental data were ac- 
curate enough to unambiguously exclude that possibility. 
It should be noted that the elastic parameters extracted 
from molecular simulation data can be severely affected 
by the short- wavelength fluctuations, especially in small 
systems. In order to avoid this problem, Fournier and 
Barbetta have carried out simulations of a discretized 
fluctuating line 'membrane' in two spatial dimensions. 
Their results not only clearly showed an effect but 
also indicated that the fluctuation tension er fluc differs 



substantially from both the internal and the frame ten- 
sion, <r lnt and <T f . The total length of the line (i.e., the 
interfacial 'area') in this study was allowed to fluctuate 
freely, only controlled by the internal tension <r int . 

Historically, the overwhelming theoretical literature on 
'tensionlcss' fluctuating fluid membranes relies on the as- 
sumption that the fluctuation tension is zero. On the 
other hand, the physically (i.e., mechanically) relevant 
tension is obviously the frame tension. If the frame ten- 
sion and the fluctuation tension are not identical, a vast 
body of work has to be revisited. Therefore, a clarifi- 
cation of the issue is clearly desirable. The purpose of 
the present work is to contribute to such a clarification. 
We will first re-analyze the theory and conclude that the 
central result, cr fluc ^ o~ t , is less rigorous than it seems. 
Then we will present numerical results which strongly 
support the hypothesis that the frame tension and the 
fluctuation tension are in fact equal for membranes with 
fixed number of lipids and fixed area. 

The starting point of the previous studies 0, [f| is 
the Monge representation for planar membranes [7| with 
fluctuating area and fixed number of lipids. The mem- 
brane is assumed to fluctuate only weakly about a 
plane, such that the local membrane position can be 
parameterized by a function h(x,y) with the constraint 
JJ dxdyh(x,y) = 0. Furthermore, the Hamiltonian is ex- 
panded in powers of h and only quadratic terms are taken 
into account. Omitting the contribution of the Gaussian 
curvature, which is constant for closed membranes with 
fixed topology, one obtains Q 
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dxdy{a {Vhf + n(Ah) 2 }, (2) 



where the integral runs of the projected area A p of the 
membrane. Since the ensemble under consideration is one 
with fixed number of lipids, er is the internal tension <7 int . 
Furthermore, the height fluctuations in this linearized 
model are clearly driven by a a , i.e., o~ nnc = a a = cr int . 

To fully define the Monge model, one needs to specify 
two additional parameters: The in-plane and perpendic- 
ular coarse-graining lengths A and A, respectively, usu- 
ally taken to be of the order of the membrane thickness. 
The parameter A acts as a small-wavelength cutoff and 
determines the number of fluctuation degrees of freedom 
TV via N = A p /A 2 . Following Refs. we will assume 
that N is proportional to the number of lipids. 

For future reference, we explicitly point out the ap- 
proximation yj\ + (Vh) 2 ps (1 + (Vh) 2 /2) entering the 
expression ([2]). The Monge model thus approximates the 
full Hamiltonian for planar interfaces up to the order 
[Vh) 2 ~{A-A P )/A P . 

Throughout the remaining paper, energies will be given 
in units of k B T. The partition function of the Monge 
model can be evaluated analytically, giving the free en- 
ergy 
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G(N, a int ,A p ) = a lat A p + i^il (lm> lnt A 2 ) - 2 (3) 



C int A p 
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Attk(N-I) . 

C i„t Ap 



The only approximation entering this result was to re- 
place the sum ^2j:f(\k\) in Fourier space k by the inte- 



gral J max dfc f(k), where the value fc max = ^JAtt^- 

is imposed by the requirement SjLso = N—l. 

Eq. ([3]) gives the intcrfacial free energy of a fluctuat- 
ing interface with fixed projected area A p and fluctuat- 
ing total area A, controlled by the internal tension cr int . 
Membranes are better described by an ensemble where 
Ap fluctuates, controlled by the frame tension er f , and A 
is fixed. In Refs. it is argued that these different 

ensembles are equivalent and the appropriate free energy 
G(N, A,a f ) can be obtained by a Legendre transform. 
Here the number of degrees of freedom TV is still fixed, 
because the number of lipids in the membrane is fixed. 
Without having to determine G explicitly, one can easily 
calculate the relation between the frame tension and the 
internal tension via a t = , giving 
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— ln(l 



int Ap 



(N-l) 



(4) 



This reproduces the result in Refs. @, 0, [H[ and demon- 
strates the existence of a residual internal tension, er int ~ 
N— 1/2 A p , at cr f = 0. Likewise, the total area is related 
to cr int and A p by 
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da int 



= A n 1 



8nn 



ln(l 



Attk(N-I) 

O 'i„t Ap 



(5) 



Inverting this relation and expressing <7; nt as a function 
of A and A p , one can rewrite the ratio a { /a iat entirely as 
a function of (A — A p )/A p , giving 



^- = 1 l -( e y-l-y) « 



(6) 



with y = 8itk(A — A p )/A p . In the last step, we have ex- 
panded t7 f /(T int in powers of the excess area {A — A p ) / A p . 

Hence the frame tension is found to differ from the in- 
ternal tension to second order in (A — A p )/A p . On the 
other hand, the Monge model approximates planar inter- 
faces only up to the first order in (V/i) 2 or (A — A p )/A p . 
Thus the results ([6]) and (j4]) arc not rigorous! Whatever 
the differences between er f , er fl n C , and er int in real mem- 
branes are, it is not possible to determine them within 
the Monge model. 

Another serious issue is the equivalence of ensembles, 
which was taken for granted in the above derivation. 
Different ensembles are only equivalent in the thermo- 
dynamic limit of infinite membranes A — > oo at fixed 
A/N. At zero tension, this limit does not exist within the 



Monge model, since tensionless membranes bend around 
on length scales larger than the persistence length and 
are no longer planar. Hence the thermodynamic limit is 
generally questionable for membranes in the floppy low- 
tension regime. 

One can also derive directly a Monge approximation 
for planar membranes in the desired (N, A, a { ) ensemble. 
The Hamiltonian would then read T-L = — a t A p + H- bondiasl 
where the second term "H b endin g includes the pure bending 
contributions, and fluctuations are subject to the con- 
straint that A is constant. Expandin g ag ain in powers of 
h up to quadratic order, one obtains |25[ 



-a, A 



dec dy{(7 f (V/i) 2 + n(Ah) 2 }. (7) 



According to this alternative Ansatz, the fluctuations are 
driven by the frame tension. This result is also only valid 
up to order (A-A p )/A p , like Eq. ©. The Monge model 
can hence be used both to predict equivalence and non- 
equivalence of (7 fl n C and cr f ! 

The question remains which is correct. To shed light on 
this issue, the present author has carried out numerical 
Monte Carlo simulations of one-dimensional 'membranes' 
in two-dimensional space. The model is similar to that 
of Fournier and Barbetta [B| , except that the membranes 
had fixed intcrfacial 'area', i.e., fixed contour length L. 
Compared to two-dimensional surfaces, one-dimensional 
membranes have the advantage that the Hclfrich Hamil- 
tonian, Eq.(Q]), can be discretized in a straightforward 
manner, without having to resort to sophisticated models 
methods like randomly triangulated sheets [2(|. More- 
over, fluctuation effects are stronger in two dimensions 
than in three dimensions, hence they can be studied more 
easily 

To set the stage, we briefly repeat the previous calcu- 
lations up to Eq. ([6]) for membranes in two-dimensional 
space. 'Membranes' are then simply lines with total 
length L and projected length L p . Otherwise, the Monge 
model is constructed in the same way as in Eq. ^ , with 
a height function h(x). The equivalent of Eq. J3]) reads 



G(N, (Ti„t, Lp) = a int L p + 

N/2 



N-l , ( \ 2 a int 2 L p 



In 
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where we have introduced the dimensionless rescaled ten- 
sion a = (7j n t L 2 /4k. We note that in two dimensions, the 
sum may not be replaced by an integral, because it 
is dominated by the contributions at small k. 

The total 'area' L and the frame tension a t can again 
be determined by taking the appropriate derivatives of 
the free energy, L = dG/da int and <r f = dG/dL p . After 



some algebra, and approximating 2fc=i(l + (^) 2 ) 
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L — Lr 



cr, 



i(xcoth(a;) — 1), one obtains 
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3(N-1) 
2L n 



(9) 
(10) 



which takes the value <r f — a iat = — 3(N— l)/2L p at a 0. 
In two dimensions, a { /a int as a function of (L — L p )/L p 
has an essential singularity at (L — L p )/L p — > 0. There- 
fore, we cannot give a rigorous expansion as in Eq. ([6]). 
but we can only derive an approximate expression for the 
regime coth(y/a) ~ 1, 
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L, 



24k- 



N-l ,L- L, n 



L, 



L,, 



(11) 



which shows that <x f /<7 int significantly deviates from one 
in two dimensions (i.e., the deviation is linear in the small 
parameter (L — L p )/L p ). 

For future reference, we also give the expression for 
the squared amplitude of fluctuations w 2 = (h 2 ) — (h) 2 , 
which was obtained by direct evaluation of the statistical 
average: 



(12) 



This completes the summary of results for the Monge 
model in two dimensions. Given the previous discus- 
sion, we would not expect them to be of much use, be- 
cause they are only valid to linear order in (L — L p )/L p 
and they have been calculated for the wrong ensemble 
((N,a int ,L p ) as opposed to (N,L,a t )). The ensembles 
can only be expected to be equivalent in the thermo- 
dynamic limit L p — > 00 at fixed N/L p , corresponding 
to a — > 00, which is a very special limit in the above 
equations. Surprisingly, it will turn out that the equa- 
tions nevertheless describe the simulation data remark- 
ably well. 

In the numerical simulations, membranes with fixed to- 
tal length L were split in N segments of fixed length L/N, 
and a frame tension af was imposed in the x-direction. 
The discretized Hamiltonian reads 



H 



07 L p + k — 2^(1 - • ej_ij 



/v 



(13) 



where the sum i runs over the segments, e, denotes the 
unit vector in the direction of the ith segment, L p is the 
projected length (L p — |f X^ e i,»)i and periodic bound- 
ary conditions were applied in the x direction. Monte 
Carlo simulations were carried out for two discretiza- 
tions, N = 40 and N = 100, and parameters k/L and 
07L ranging from k/L = 1 to k/L = 1000 and OfL = 
and GfL = 0.05 to 07L = 10000. Typical run lengths 
were 0.5 — 1.0 • 10 9 Monte Carlo steps. Since the re- 
sults are compared with the predictions of the Monge 
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FIG. 1: Difference between internal tension a int and frame ten- 
sion crfvs. frame tension a t . Inset: Raw data for discretization 
N = 100 and selected values of the stiffness k. Main frame: 
Rescaled data for two discretizations iV and all values of k 
(k/L = 1, 2, 3, 5, 10, 15, 20, 30, 50, 100, 200, 300, 500, 1000), vs. 
rescaled tension a = a/(L 2 /4«), compared with the theoreti- 
cal master curve f(x) — v / a;coth(Y / s) — 1 from Eq. (|10[) 



model, parameter combinations (k,<t/) were disregarded 
for which more than 1 / 1000 of segments had angles larger 
than 45° with the x-axis. This applied to k/L < 1 at 
frame tensions afL < 30. From the simulation data, 
one can extract in a straightforward manner the sta- 
tistically averaged projected length, (L p ), the squared 
amplitude of fluctuations, (w 2 ), and the internal ten- 
sion cr int = -dG/dL = -(&H/dL) + 2(JV-lVi, where 
G(N, L, <7 f ) is the free energy of the system [271 ] 

Fig. [U summarizes the results for the internal tension. 
The bare data in the inset clearly show that the inter- 
nal tension deviates from the applied frame tension, as 
predicted by Eq. (fT0|) . The value of the residual tension 
at cr f — > however disagrees with the theory. Whereas 
Eq. pTJ)) predicts a int L\a f =o ~ 3(7V— 1)/2, the actual data 
rather suggest <7 int L\cr f =o = 5(JV— 1)/2. The additional 
term (N—l)/2 can be attributed to the translational free- 
dom of the segments along the x-direction. 

Next we test whether the data for (cr f — a int )L p + 
3(iV— 1)/2 can be scaled in a way that they collapse onto 
a single curve for all values of k, er f , and N. Eq. (|10l) sug- 
gests to try the scaling variable a = a int L 2 /4K. It turns 
out that the data do not scale at all with this Ansatz. The 
scaling is much better with the scaling variable a t L 2 /AK, 
but still not perfect (data not shown). If one however re- 
places L p by L everywhere and uses the scaling variable 
a = tjf L 2 /4k, the data collapse almost perfectly onto 
the theoretical master curve suggested by Eq. (fTO)) for all 
values of k and af and both discretizations N . Only for 
large values of a does one observe deviations, which can 
be attributed to discretization effects since they become 
smaller for larger N (Fig. [I] main frame). 

These first findings already suggest that the relevant 
quantity driving the fluctuations is the frame tension, 




o f L 2 /4k 

FIG. 2: Scaling plots of the squared amplitude of fluctuations 
w 2 vs. rescaled tension a = <jf{L 2 /An) compared with the 
theoretical master curve f(x) = ||(f + 1 — \fi~ coth(y / x)) 
(Eq. (|12[) ~). Inset shows an alternative scaling plot where 07 
in a is replaced by a' = af + 1/L. 
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FIG. 3: Scaling plot of the excess length {(L — L p )), vs. 
rescaled tension a = a f(L 2 / '4k), compared with the theo- 
retical master curve, f(x) = ^(y / xcoth(y / x) — 1) (Eq. (©). 
Inset compares the scaling with shifted and unshifted a' in 
the most critical parameter region. 



not the internal tension. Furthermore, they also suggest 
that the relevant length scale in the system is the total 
length L, not the projected length L p . With the cor- 
responding replacements, the results are in surprisingly 
good agreement with the predictions of the Mongc model, 
despite the fact that the latter have been derived for a 
different ensemble with different scaling variables. We 
will now proceed to investigate two quantities that probe 
the fluctuations directly, i.e., the squared amplitude of 
fluctuations w 2 and the excess length (L — L p ). These 
quantities characterize the fluctuations in an integrated 
way. 

Fig. [2] shows scaling plots for the squared amplitude 
of fluctuations. The data are rescaled as suggested by 
Eq. (|T2"j) . with L p replaced by L and plotted against the 
rescaled frame tension a. The scaling is quite good, but 
not perfect. It can be improved by shifting a c by the small 
amount 5a = 1/L (see inset of Fig. [2]) . The data then col- 
lapse very nicely onto the theoretical curve predicted by 
(fT2]l . No data collapse is obtained when plotting against 
the rescaled internal tension cr int (not shown) . Thus the 
results for the squared amplitude of fluctuations indicate 
that the fluctuations are driven by the frame tension cr f , 
possibly slightly shifted, and not by the internal tension 
<r int . The results for the excess length (L — L p ) support 
this conclusion (Fig. ([3])). In that case, the best scaling 
is achieved directly with a = <jf (L 2 /4k), without addi- 
tional shift. The data again collapse excellently onto the 
theoretical master curve suggested by Eq. (j9|). 

In sum, the simulation data show clearly that the 
membrane fluctuations are driven by the frame tension 
<7 f . These results raise the question why Barbetta and 
Fournier in Ref. [j| obtained a different result from sim- 
ulations of essentially the same model. The answer is 
that these simulations were carried out in a different en- 
semble where both the contour length and the projected 



length were allowed to fluctuate, i.e., the (TV, <7 int , cr f ) en- 
semble as opposed to the (TV, L, <j f ) ensemble considered 
in the present work. As pointed out earlier, the ensem- 
bles are not equivalent. To further analyze the issue, 
additional simulations were carried out for membranes in 
the (TV, L, L p ) ensemble and in the (TV, cr int , L p ) ensemble. 
The purpose of these simulations was to find out whether 
the fluctuations in other ensembles might be driven by 
(Ti nt instead of a { . 

Fig. 2] shows the squared amplitude of fluctuations w 2 
vs. the excess length (L — L p ), in a scaling plot which is 
independent of the actual choice of the rescaled tension. 
The data are taken from simulations in the (TV, L, a t ) en- 
semble and in the (N,L,L P ) ensemble. In the (N,L,a t ) 
ensemble, they collapse onto a single curve which agrees 
nicely with the theoretical prediction from Eqs. © com- 
bined with (fT2]l . The data from the (N,L,L p ) ensemble 
also collapse, but the scaling curve is different and not 
consistent with Eqs. © and |(T2J) [H. 

If L is allowed to fluctuate, the situation is further 
complicated by the fact that in one dimension, the to- 
tal bending energy of bubbles decreases with increasing 
bubble size. As a consequence, the simulation data in 
the (TV, cr int , L p ) ensemble feature a transition between a 
flat state and an inflated bubble state at bending stiff- 
ness k/L ~ 3 and low tensions (data not shown). The 
simulations of Fournier and Barbetta were carried out at 
k/L = 2.5 (in yet another ensemble), hence the vicinity 
of a similar transition possibly accounts for some of the 
intriguing phenomena reported in Ref. 0. 

In the introduction, we have raised the question 
whether theories and simulations / experiments really con- 
sider the same object if they refer to 'tensionless' mem- 
branes? Based on the simulation results presented above, 
we conclude that the answer is most likely "y es " • At least 
in one dimensional membranes with fixed 'area', the flue- 
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FIG. 4: Rescaled squared amplitude of fluctuations vs. excess 
length. Circles and stars correspond to data for fluctuating 
projected length L p and discretizations N = 40 and N = 100, 
respectively. Diamonds correspond to data from simulations 
with fixed projected length L p . The solid line shows the the- 
oretical prediction from Eq. <(9j combined with Eq. (|12[) . The 
dashed line is a guide for the eye. 

tuations are driven by the imposed lateral stress, and the 
behavior of stress-free (tensionless) membranes is consis- 
tent with that of the Mongc model at tension zero. These 
results will presumably also hold for two-dimensional 
membranes. We have presented an argument in Eq. flTJ) 
that rationalizes the prediction a anc = o~ t in the physically 
relevant (N, A, er f ) ensemble for arbitrary dimensions. We 
should note, however, that our numerical simulations in- 
dicate that the relevant length scale is the total 'area' 
rather than the projected 'area'. This presumably renor- 
malizes the actual 'g 2 coefficient', er fluc = a t (L p /L) 2 in 
two dimensions (29j . or o~ auc = a t (A p /A) in d dimen- 
sions. A detailed analysis of the full fluctuation spectra 
as a function of the wavevector q will be presented else- 
where. 

The present author hopes that this paper will stim- 
ulate discussions and further work on these intriguing 
and important issues. For example, it is not clear why 
the Monge model provides such a good prediction of the 
scaling functions for the various quantities, even though 
it is derived for the wrong ensemble with different scal- 
ing variables. Furthermore, mesoscopic simulations of 
two-dimensional fluctuating fluid sheets, e.g., using ran- 
domly triangulated surfaces [2fjj or Fourier techniques 
[I(J , would clearly be desirable. 

The author thanks J. Neder, J.-B. Fournier and M. 
Weikl for inspiring discussions. Partial support from the 
German Science Foundation (CRC 625) is gratefully ac- 
knowledged. 



New York, 2002). 
[2] O. Farago, P. Pincus, Eur. Phys. J. E 11, 399 (2003). 
[3] O. Farago, P. Pincus, J. Chem. Phys. 120, 2934 (2004). 
[4] A. Imparato, J. Chem. Phys. 124, 154714 (2006). 
[5] J.-B. Fournier, C. Barbetta, Phys. Rev. Lett. 100, 

078103 (2008). 
[6] J. Stecki, J. Phys. Chem. B 112, 4246 (2008). 
[7] S. Safran, Statistical Thermodynamics of Surfaces, Inter- 
faces, and Membranes (Westview Press, Boulder, 2003). 
[8] W. Helfrich, Z. Naturforschung 28c, 693 (1973). 
[9] J. S. Rowlinson, B. Widom, Molecular Theory of Capil- 
larity (Clarendon, Oxford, 1982). 
[10] The fluctuations of the height h(x,y) are expanded as 

(\hq\ 2 ) = k B T/{o ail , : q 2 + 0(q 4 )) in Fourier space q. 
[11] W. Cai, T.C. Lubensky, P. Nelson, T. Powers, J. Phys. 

II (France) 4, 931 (1994). 
[12] E. Evans, W. Rawicz, Phys. Rev. Lett. 64, 2094 (1990). 
[13] P. Girard, J. Prost, P. Bassereau, Phys. Rev. Lett. 94, 
088102 (2005). 

[14] K. Sengupta, L. Limozin, Phys. Rev. Lett. 104, 088101 
(2010). 

[15] R. Goetz, G. Gompper, R. Lipowsky, Phys. Rev. Lett. 

82, 221 (1999). 
[16] E. Lindahl, O. Edholm, Biophys. J. 79, 428 (2000). 
[17] S.J. Marrink, A.E. Mark, J. Phys. Chem. B 105, 6122 

(2001). 

[18] Z.J. Wang, D. Frenkel, J. Chem. Phys. 122, 234711 

(2005) . 

[19] W.K. den Otter, J. Chem. Phys. 123, 214906 (2005). 
[20] G. Brannigan, F.L.H. Brown, Biophys. J. 90, 1501 

(2006) . 

[21] B. West, F.L.H. Brown, F. Schmid, Biophys. J. 96, 101 
(2009). 

[22] C. Loison, M. Mareschal, K. Kremer, F. Schmid, J. 

Chem. Phys. 119, 13138 (2003). 
[23] J. Stecki, J. Chem. Phys. 125, 154902 (2006). 
[24] J. Neder, PNielaba, B. West, F. Schmid, J. Chem. Phys. 

132, 115101 (2010). 
[25] Note that the base A p of the integral f A fluctuates. To 

restore a standard Monge model with fixed base, one may, 

e.g., replace f A by / ( ^ by rescaling x ->• x y/(Ap)/A P , 

y — > y \J (A p ) /A p in a self-consistent manner. This will 
only introduce another correction of order ((A—A p ) /A p ) 2 
in the Hamiltonian H' Monga . 
[26] G. Gompper, D.M. Kroll, J. Phys.: Cond. Matter 9, 8795 
(1997). 

[27] This is analogous to the virial equation for the pressure in 
simple (NVT) systems. To derive it one rewrites the par- 
tition function Z = Ud =2 [f dri<5(|r; ~ ~ L/N)]e~ n 
with the joint positions r*i (fixing fb = 0) in terms 
of rescaled variable fj = fj L/N and then calculates 
-dG/dL = d{\nZ)/dL. 

[28] The agreement with the theory does not become better 
if one rescales w 2 and (L — L p ) with L p instead of L. 

[29] For a given wave number n, the height fluctuations are 
(h n h-n) oc 1/(t n 2 + 0(n 4 )) with r n 2 = cr t (2-Kn/L) 2 = 
o auc (2im/L p ) 2 . 

[30] A. Troster, Comp. Phys. Comm. 179, 30 (2008). 



[1] B. Alberts, A. Johnson, J. Lewis, M. Raff, K. Roberts, P. 
Walter, Molecular Biology of the Cell, 4th ed. (Garland, 



